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Abstract 

o 

Contribution of matter fields to the Gell-Mann-Low function for N=l supersym- 
metric Yang-Mills theory, regularized by higher covariant derivatives, is obtained 
using Schwinger-Dyson equations and Slavnov-Tailor identities. A possible devia- 
tion of the result from the corresponding contribution in the exact Novikov, Shifman, 
Vainshtein and Zakharov (3- function is discussed. 

CD . 

1 Introduction. 

> 
X 

It is well known, that supersymmetry essentially improves the ultraviolet behavior of 
a theory. For example, even in theories with unextended supersymmetry, it is possible to 
propose the form of the /3-function exactly to all orders of the perturbation theory. This 
proposal was first made in Ref. [T] as a result of investigating instanton contributions 
structure. For the N — 1 supersymmetric Yang-Mills theory with matter fields, this 
/3-function, called the exact Novikov-Shifman-Vainshtein-Zakharov (NSVZ) /3-function) 



is 



/3(a) 



a 2 



3C 2 - 2C(R) (l - 7 (a) 



2tt(1 - C 2 a/2n) ' v ' 

where 7(a) is the anomalous dimension of the matter superfield. We note that here 
we present the result for a theory with the Dirac fermions, and two chiral superfields 
correspond to each of them. R denotes a representation for one of the matter superfields, 
and C(R) is given by 

tr (T a T b ) = C{R) 5 ab . (2) 
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In order to obtain such /3-function it is most convenient to use the higher derivative 
regularization E]- Then the renormalization of the operator W a C ab Wb is exhausted 
at the one-loop, while the Gell-Mann-Low function coincides with the exact NSVZ (3- 
function. This was confirmed by the explicit calculations in two-loop [HIE], three-loop pj, 
and partially four-loop approximations [7j. (We note that the Gell-Mann-Low function 
does not depend on a particular choice of the renormalization scheme due to its definition.) 

Investigation of the N — 1 supersymmetric electrodynamics was made in Ref. |H] 
exactly to all orders of the perturbation theory. According to that paper, it is possible 
to obtain the exact /3-function using Schwinger-Dyson equations and Ward identities. 
However, there is an additional contribution besides the ordinary NSVZ result. Explicit 
calculations show that it is always equal to in the lowest orders. This allows proposing 
existence of a new identity for some Green functions. A proof of this identity exactly to 
all orders is rather nontrivial. The matter is that it is not reduced to the Ward identity, 
while its proving by summation of diagrams |Hj is rather complicated. We note that in 
Ref. this proof was made only for a restricted class of diagrams in the massless case. 
Nevertheless, it seems that the new identity is also valid for the other diagrams, that 
follows, for example, from the explicit four-loop calculations [Zj. 

So, obtaining the exact /3-function is a rather nontrivial problem even in the elec- 
trodynamics. In the case of the Yang-Mills theory the exact /3-function agrees with the 
results of explicit calculations in the two-loop approximation if the dimensional reduction 
is used for regularization. However, there are disagreements between the predictions, 
following from the exact /3-function, and explicit calculations, made in MS-scheme for the 
/5-function, defined as a derivative of the renormalized coupling constant, in the three-loop 
approximation [El EH [TJ] . It was shown [TT] that these disagreements can be eliminated 
by a special choice of renormalization scheme, the possibility of such a choice being highly 
nontrivial |1 3| . 

In this paper we attempt to investigate contribution of the matter superfields to the 
Gell-Mann-Low function in non-Abelian supersymmetric gauge theories exactly to all 
orders of the perturbation theory using the Schwinger-Dyson equations and Slavnov- 
Taylor identities. The Gell-Mann-Low function is defined by asymptotic of the two-point 
Green function in the limit of large momentums. And so, the obtained results actually 
correspond to the limit of massless matter fields. (Nevertheless, in order to elucidate 
calculating the contribution of Pauli-Villars fields, we will consider the case of an arbitrary 
mass throughout the paper. Taking the limit m — > will be made in the end for obtaining 
the final result.) 

This paper is organized as follows. 

In Sec. El we recall basic information about the N — 1 supersymmetric Yang-Mills 
theory and its higher derivatives regularization. The Schwinger-Dyson equations for the 
considered theory are obtained in Sec. 03 In that section we present the solution of 
the Slavnov-Taylor identities, which determines vertex functions in the Schwinger-Dyson 
equations. After substitution of these vertexes the contribution to the two-point Green 
function of the gauge field is obtained in the form of an integral. Calculation of this 
integral is made in Sec. |H A brief discussion of the results is given in the conclusion. 
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2 N = 1 supersymmetric Yang-Mills theory, back- 
ground field method and higher derivative regular- 
ization 



The N = 1 supersymmetric Yang-Mills theory in the superspace is described by the 
action 

S = ^Retr J d 4 xd 2 9W a C ab W b + d 4 xd 4 9 (Ve 2y + 4> + e~ 2Vt ^j + 
+-m [ d A xd 2 6 4> t <f) + -m [ d 4 xd 2 64> + <P*. (3) 



2 J T T 2 

Here and are chiral matter superfields, and V is a real scalar superfield, which contains 
the gauge field as a component. The superfield W a is a supersymmetric analogue of 
the gauge field stress tensor. It is defined by 



(4) 



where 



D = Z.-vfOd li (5) 

is a supersymmetric covariant derivative. In our notation, the gauge superfield V is 
expanded over the generators of the gauge group T a as V = e V a T a , where e is a coupling 
constant. 

Action is invariant under gauge transformations 

- e iA 0; ^ e" 4At 0; e 2y ^ e iA+ e 2y e - jA . (6) 

Such a transformation law means that if the field <p is in the representation R of the gauge 
group G, the field is in the representation R, conjugated to R. 

For quantization of this model it is convenient to use the background field method. 
The matter is that the background field method allows calculating the effective action 
without manifest breaking of the gauge invariance. In the supersymmetric case it can 
formulated as follows [HJUn]: Let us make a substitution 

e w ^ e 2V = e^+e^e", (7) 

in action |JBJ), where O is a background scalar superfield. Expression for V is a complicated 
nonlinear function of V, f2, and Q + . We do not interested in explicit form of this function: 

v' = v\y,n]. (8) 

(For brevity of notation we will not explicitly write the dependence on fl + here and 
below.) 

The obtained theory will be invariant under background gauge transformations 
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e iK e%~ A ; 



V -> e 4A 1/e 
e- A+ e" + e-^, 



(9) 



where K is a real superfield, and A is a chiral superfield. 
Let us construct the chiral covariant derivatives 



D 



= e 



(10) 



Acting on some field X, which is transformed asX-> e lK X, these covariant derivatives 
are transformed in the same way. It is also possible to define a covariant derivative with 
the Lorentz index 



D, 



(crr b {D a ,D b y, 



which will have the same property. It is easy to see that after substitution (£[J) action Q 
will be 



S= ArtrRe ! d 4 xd 2 9W a W a ^rtrRe / d 4 xd 4 9 

Ae 2 / 128e 2 



\Q(e- 2V D a e 2V )w a 



+ {e- 2V D a e 2V )D 2 e- 2V D a e 2V 



1 / J^/'tflo+r^V-' <" -<r< ~ u <-- X 



xe 



-Si 



1 x V^rn I d 4 xd 2 9(j) t (f)+^m / <I l .rJ 2 f) t 



(12) 



where 



and the notation 



±e a D(l - 75 )^(e-"e-" + (l + 7 5)^e" + e n ) e - !J 



(13) 



r> 2 



if)(l+ T5 )r>; J D 2 = iD(l- 7 5)A 



L ^(l + 75)' 
^(l-7s) 



D, 



D a 



1 

2 1 



(14) 



is used. Action of the covariant derivatives on the field V in the adjoint representation 
is defined by the standard way. 

It is convenient to choose a regularization and gauge fixing so that invariance © 
will be unbroken. For example, the higher covariant derivative regularization can be 
introduced by the substitution of action (|T2"1) for 



4 



S A = -I-trRe [ d A xd 2 9W a W a ^rtrRe / d A xd 4 6 

Ae 2 / 128^ 



+ e 



-* v D a e* v ) D 2 {e~' v D a (l + ) e 2V 



+ ^ I d A xd 4 e 



(15) 



This action is evidently invariant under background gauge transformations (jHJ). 

We note that the described way of regularization is a bit different from the method, 
proposed in Ref. ^H]. The difference is in the form of the term, which contains higher 
derivatives. Using of the method, proposed here, simplifies calculations in a certain degree, 
while all particular features of higher derivative regularization are the same in the both 
cases. So, although explicit calculations by the method, proposed in Ref. [IB] in two 
and more loops were not yet made, it is possible to propose that their results will not 
essentially differ from the ones, obtained below in this paper. 

The gauge can be also fixed by the invariant way, for example, by adding the terms 



S " = "6i? / Me(vD^{l + { -^-)v + VD^ (l + ^)f). (16) 
In this case terms quadratic in the superfield V will have the simplest form: 

1 r / (D 2 ) n \ 

-trRejd^d A eVDl{l + [ -^L)v. (17) 

Certainly, it is also necessary to add the corresponding action S g h for the Faddeev-Popov 
ghosts. We will not write here explicit form of this action, because it is not essential for 
calculating a contribution of the matter superfield loops. 

Proposed way of the regularization and gauge fixing preserves both invariance under 
the supersymmetry transformations and the invariance under transformations (JHJ). As 
a consequence, the effective action, calculated by the background field method, will be 
invariant under both supersymmetry and gauge transformations. 

Let us construct the generating functional as follows: 

Z[J, Q, j] = J Dfi exp {iSl + iSgf + iS gh + iS s + iS^ + 

+i J d 4 x d 4 e (j + j[n]) (y'[v, n] - v) } , (is) 

where the superfield V is given by 

e 2V = e H+ e n (1Q) 

and J[fl] is a so far undefined functional. A reason of its introducing will be clear later. 
The functional integration measure is written as 
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Dfi = DV D<pD(j). 
S A denotes renormalized and regularized action 



(20) 



S r A = ^trRe f d 4 xd 2 6W a W a - T^rtrRe / r/'.rr/'fl 



4e 2 



128e 2 



+ (e- 2V D a e 2V )D 2 (e- 2V D a (l + 



D 



A 2 ™ 



/" \ e 2V 



w(e- 2V D a e 2V )w a + 

^Z(e,A//i) J d 4 xd 4 9 (Ve" + x 

/>*, (21) 



/ d 4 x d 2 9 0* + / r/\rr/ 2 f9 



where eo is a bare coupling constant, and Z is a renormalization constant for the wave 
function of the matter superfields. S g f denotes gauge fixing terms (JT6l) (the coupling 
constant e in them is also replaced by eo), and S g h is a corresponding action for the 
Faddeev-Popov ghosts (and also for the Nielsen-Kallosh ghosts). S$ denotes the terms 
with sources for chiral superfields, which in the extended form are written as 

s s = J d 4 x d 2 e (> + ~f 0) + J d A x d 2 e (j + 0* + j+0*) . (22) 

Moreover, in generating functional (JT%1) we introduce the supplementary sources 



S(h = \J d 4 x d 4 9 (0 o + e n+ e 2y e n + + e"V y e"0 o + 



+0ne 



-O* -2V* -si* 



+ 



,-n* -2V* -n* 



A) I, 



(23) 



where 0o, 0q , 0o and 0g are arbitrary scalar superfields. In principle, it is not necessary 
to introduce the term in the generating functional, but the presence of the parameters 
0o is highly desirable for investigating the Schwinger-Dyson equations. 

In order to understand how generating functional (JT8|) is related with the ordinary 
effective action, we perform the substitution V — > V . Then we obtain 



z[j,sij] = ex P { -i J d 4 xd 4 e (j + j[n]Jv}z j + j[n],n,j 



(24) 



where 



Z [J, n,j] = j Dfi exp ^iS r A + iS g f + iS gh + iS s + iSfo + i j d 4 xd 4 9JV^. (25) 

If the dependence of S g f, S g h, and S<f, on the arguments V, fl, and fl + were factorized 
into the dependence on the variable V, Z Q would not depend on Q and fl + and would 
coincide with the ordinary generating functional. This really takes place for action Q and 
Sfa. However, in the term with the higher derivatives and in the gauge fixing terms such 
factorization does not occur. Therefore, Z actually differs from the ordinary generating 
functional. 



6 



Using the functional Z[J, it is possible to construct the generating functional for 
the connected Green functions 



w[j,n,j} = - i \nz[j,n, 3 ] = 

d A xd A e(j + j[o])v + w j + j[n],n,j 



Also it is possible to construct the corresponding effective action 



(26) 



r[v,n,<f>) = - J d 4 xd 4 9(jv + j[n]v) + w j + j[n),nj 

-J d 4 xd 4 6JV- J d 4 xd 2 e( K i t + - J d 4 xd 2 e(j + <p* (27) 
where the sources should be expressed in terms of fields using the equation 



v = jjW[j,n,j] = -v + jjW 



jrW[j,n,j] = j,w j + j[n],n,j 



j + j[n},n l3 

e.t.c. 



Substituting these expression into Eq. ([27]) . we write the effective action as 



(28) 



r[v,n,. 



Wn 



j + j[n],n,j 



d*xd*Q (j\ft]V + j— w 

V o J 

similar terms with d> + , < 
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d 4 xd 2 6<j)—W J+J[Q],n,j 



j,j + j[n},n,j 

< and d> 



(29) 



Let us now set V = 0, so that 



j+j[G\,n,j 



(30) 



Moreover, we take into account that the invariance under background gauge transfor- 
mations (JOJ essentially restricts the form of the effective action. If the quantum field V in 
the effective action is set to 0, the superfield K will affect only the gauge transformation 
law of the fields f2 and f2 + , and the only invariant combination is expression (fT9|) . (It 
is invariant in a sense, that the corresponding transformation law does not contain the 
superfield K.) This means that in the final expression for the effective action we can set 



V. 



(31) 



In this case the effective action is 



r[o,v, 



Wo 
5 



J+J[V],V,j - / d 4 xd 4 6 ( J + J[Y))-^W 



- / d*xd z e<j)— W 



J + J[V),V,j 



/ 5 J 

similar terms with 



J + J[V],V,j 
h, and 4 



(32) 



7 



Note, that this expression does not depend on form of the functional J[ft]. In particular 
it can be chosen to cancel terms linear in the field V in Eq. (fTHl) . Such a choice will be 
very convenient below. 

If the gauge fixing terms and the terms with higher derivatives depended only on 
V, expression |32l) would coincide with the ordinary effective action. However, as we 
already mentioned above, the dependence on V, fl, and f2 + is not factorized into the 
dependence on V' in the proposed method of renormalization and gauge fixing. According 
to Ref. [13 HH] the invariant charge (and, therefore, the Gell-Mann-Low function) is 
gauge independent, and the dependence of the effective action on gauge can be eliminated 
by renormalization of the wave functions of the gauge field, ghosts, and matter fields. 
Therefore, for calculating the Gell-Mann-Low function we may use the background gauge 
described above. We note that if this gauge is used, the renormalization constant of the 
gauge field is 1 due to the invariance of the action under transformations ©. 

Nevertheless, generating functional (fT8|) is not yet completely constructed. The mat- 
ter is that adding the term with higher derivatives does not remove divergences from 
one-loop diagrams. To regularize them, it is necessary to insert the Pauli-Villars deter- 
minants in the generating functional [HU- Because in this paper we are interested only in 
contributions of matter superfields, we construct these determinants only for them. (Also 
it is necessary to introduce the Pauli-Villars fields for the gauge field and ghosts.) Here 
we will at once use condition (|3T| . So, we insert in the generating functional the factor 

H^detPV^V^^y, (33) 

i 

in which the Pauli-Villars determinants are defined by 

(detPy(V,M)) l = J D$D§ exp (iS PV ^. (34) 
If condition (|3TJ) is used, the action for the Pauli-Villars fields will be 

S PV = Z(e, A/p)i j d A xd A 6 ($ + e v e 2V/ e v $ + $+ e - v V 2V V vt <!>) + 
+^ J d 4 xd 2 6M&$ + ^ J d 4 xd 2 6M$ + $*. (35) 
The coefficients q in Eq. (j33|) satisfy conditions 

J> = 1; J>M 2 = 0. (36) 

i i 

Below, we assume that M, = OjA, where Oj are some constants. Inserting the Pauli- 
Villars determinants allows cancelling the remaining divergences in all one-loop diagrams, 
including diagrams containing counterterm insertions. 
Making the transformation 

<!>^<5>/VZ; I -> Q/y/z (37) 

in the generating functional, we obtain that the dependence on the renormalization con- 
stant Z can be easily found from the generating functional, which does not contain the 
constant Z. If <p and j are 0, this is made by the substitution 
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m^m/Z- Mi/Z. (38) 

Therefore, we will make calculations with Z — 1, and dependence on the renormalization 
constant will be restored in the final result. 

3 Schwinger-Dyson equations and Slavnov-Taylor 
identities. 



In order to construct the Schwinger-Dyson equations for the considered theory it 
is necessary to split the action into three parts: the action for the background field, 
the kinetic term for quantum fields, which does not contain the background field, and 
interaction, in which the other terms are included: 

S = S(V) + S 2 (V, 0) + Sj(V, V, 0). (39) 

(Earlier we saw that the terms of the first order in the superfield V, which were obtained 
from the expansion of the classical action, can be omitted.) So, generating functional (fT8|) 
can be written as 



Z[J,j, V] = J dfj, exp (iS[V) + iS 2 [V, <p] + iS T [V, V, 0] + iS s + i J d s x JV^j = (40) 



exp (iS[V]+iS_ 



15 15 15 



i 5 J i 53 % 5j 
x J dfx exp (iS 2 [V, 4>]+iS s + i J d 8 x JV + i J d 8 x JV 



j=o 



where J d^x = / d^xd 4 6 x . Let us differentiate this expression with respect to the back- 
ground field 



d Z[J, j, V] = Z + exp ( iSi 



5Y X 



5V 



15_ 15_ 15_ 

i 5 J 1 i 5 J ' i 5j 



iJ x x 



x J dfx exp (iS 2 [V, <p] + iS s + i J d s x JV + i J d 8 x jvj 



(41) 



j=o 



Moving the current J x to the left and dividing the result to Z, we obtain 

5 



SV. 



% 5 J i 5j\ 



(42) 



Because the background field V is a parameter of the effective action, these equality can 
be equivalently written as 



sv _ ssyv] + ^ 



5V 



5V X 



5Vn 



(43) 
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where the angular brackets denote taking an expectation value by the ordinary functional 
integration. Certainly, it is necessary to set the field V (an argument of the effective 
action) to in the final result in accordance with Eq. (|30|) and the corresponding discus- 
sion. 

Let us find a contribution to this expression given by the matter superfields. The 
corresponding interaction terms are 

Sj = ^J d 8 x <p + (e v e 2V e v - lj <j> + similar terms with 0. (44) 

Differentiating Si with respect to the background field, we obtain that the corresponding 
contribution to the effective action is written in the form 



— — = terms without matter + V a— — ( In det Pv(v, V, Mj) ) + 



+ --^— J d 8 x (^(p^e Vx e 2Vx e Vx (f) x + similar terms with 0y. (45) 

We are interested in the two-point Green function of the gauge field, corresponding to 
the expansion of the effective action in powers of the background field up to the second 
order terms. To calculate such expansion, it is convenient to use the following equality 



V 2V V 

e e e 



v=o 



V i Ye 2V + e 2V v 



(y e 2V + e 2V V)V 



e / d 8 xd 8 yV b y V a x 



v 2V v 

e e e 



v=o ™ y 



V 2V V 

e e e 



rpa 



V=0 



(46) 



where e in the beginning of the second string is a coupling constant. Using the form of 
the action for supplementary sources (|23|) . we easily obtain 



1 

-tr 

l 



5 2 T S 2 T 
+ 



(47) 



(In the last equation the matrix notation is used for the brevity.) The derivatives with 
respect to the sources must be expressed in terms of fields as follows: 



d 8 z 
5 2 Y 



5 2 T 



-i 



D 2 5 



+ 



5 2 r 



D' 1 



D 2 S 



8d 2 



— + 



5 2 r 



-1 



D 2 6 



8d 2 



8d 2 



+ 



+ 



5 2 T 



5V r , 



-i 



5 



(48) 



Therefore, using Eqs. (JUJ), (|4*H|l. and taking into account similar terms with the fields 
(f), the corresponding contribution to the two-point Green function of the matter superfield 
can be written as 
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s 2 r 



+ e 



5 



tr 



5 2 r 



- ( T 



£ 2 r 

^5* 



h.c. 



(49) 



i 5j+ 54b - ^ ^rox 

where dots denote contributions of the gauge fields, ghosts, and also all possible Pauli- 
Villars fields. (In this expression we omit symmetrization with respect to the indexes a 
and b, because, as we will see below, it will be symmetric automatically.) We note that 
the calculation of the Pauli-Villars fields contributions are made completely similar to 
the calculation of ordinary fields contributions (Hj, and details of this calculation are not 
presented here. The result will be given below. 

In order to calculate expressions in Eq. ijlHj). it is necessary to use Eq. ijlHJl and 
manifestly perform the differentiation. Then the result can be graphically presented as a 
sum of two effective diagrams 



Ar 



(2) 



-o 




The double lines correspond to the effective propagators, which are written as 



5 2 T 



-l 



GDlDl 



A(d 2 G 2 + m?.pf xy} 



5 2 T 



-1 



d 2 G 



mJDl 



! J 2 



xyi 



(50) 



(51) 



depending on the chirality of the ends. The functions G and J are determined by the 
two-point Green functions of the matter superfield as 



5 2 T 



D 2 D 2 
x X -G{d 2 )5l- 



5 2 T 



16 



^mJ(d 2 )5 8 
4 



xi)- 



(52) 



where 5* y = 5 (x — y)5 (6 X — 9 y ), and the subscripts denote points where the expressions 
are evaluated. 

The vertex functions can be obtained from the Slavnov-Taylor identities similar to 
the case of the electrodynamics. Certainly, it is necessary to take into account that we 
consider the vertex functions, which have an external line of the background field. As 
we already mentioned above, the effective action is invariant under the background gauge 
transformations. It is easy to see that this invariance can be expressed by the equality 



5T 5T D 2 ^ , D 2 5T 5T ^ 



5V a 

V 



+similar terms with 



5cj) y 8d 2 
and 0o 



>0y 



(53) 



where 



A, 



A 



x90x, 



^Ox — ^Ox-^-x > 



*V ai = -i(A« + A+)+0(V). 



(54) 
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Here A is an arbitrary chiral superfield, and all other terms in <5V are proportional at 
least to the first degree of the background field. (For the fields <f> e.t.c. the transformation 
laws can be also easily written.) Let us differentiate Eq. ijHHj) with respect to Ay, if)+, (fi x 
or with respect to A" ip z , <f) x . As a result we obtain the Slavnov-Taylor identities 



where 



i tf 3 r $ 2 t / $ 2 t \ 

= ±(D 2 + D 2 ) — — — T a D 2 S 8 -D 2 5 8 T a V 

= ^(D 2 + D 2 ) £J—T*tiW + D 2 (g* M T*-JJ—) . (55) 



sikr- 1 ^^ (66) 

They differ from the similar identities for the electrodynamics only in the presence of 
gauge group generators. Therefore, the solution will be also similar 



5 3 T 



5V*84$JL 



p=0 



(57) 



5 3 T 



p=0 



-D 2 D 2 6 8 6 8 

y y xy yz 



+ 



mfj'(q 2 ) J(q 2 )-1 



16 



D 2 D 2 5 8 D2yD H 8 + D 2 D 2 5 8 q 2 5- 



y y xy 



16 



The primes denote derivatives with respect to q 2 , 

1 



n V2 



D a D 2 D n 



1 -D a D 2 D n 



rpa 

(58) 
(59) 



g d 2 " 8 Q2- 

is a supersymmetric transverse projection operator, notation for the derivatives are similar 
to lfT4j) . and the functions F, f and h can not be determined from the Slavnov-Taylor 
identities. 

Functions (|H7I) and allow finding ordinary Green functions by the identities 



Di 



S 3 T 



6 3 T 



Di 



5 3 T 



5 3 T 



2 SV'SfaSfa 5V26<l> v 6<i>+' 2 5\-5<p y 5^ z SW%S<f> v 8^ z ^ 

We note that the supplementary sources in Eq. were specially introduced in order 
that such identities take place. A proof of these relations can be made completely similar 
to the case of the electrodynamics, considered in [B] in details. Using Eqs. (|60"|) we find 
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5 3 r 



5 3 T 



p=0 



Fn 1/2y [Dl5» y D*8» z )F( q 



+^G\ q 2 )DYi 5 D y [Dy xy Dy yz ) + l -Di5i y Dy yz G{q 2 ) 



p=0 



-— J'(g 2 ) 
32 w ' 



D 2 5 8 D 2 D 2 5 8 -D 2 D 2 5 S D 2 5 s 



rpa 



T a ; (61) 
(62) 



Obtained expressions for the vertex functions allow calculating all expressions, entering 
Eq. ijlHJ) in the limit of zero external momentum. For this purpose it is convenient to make 
the following observation: The first diagram in Eq. (|50|) can be equivalently presented in 
the form 




+e 



8 



5V b 

y 



tr 




T a mD l 



16d 2 V 50+50+ ; z=z 16d 2 
The single line in the second diagram corresponds to 

16d 2 12 



5 2 T 



-i 



+ h.c. 



(63) 



(64) 



if fields at the ends 1 and 2 have opposite chirality, and 0, if chirality of the ends is 
the same. Eq. (|63l can be easily verified using Feynman rules and performing the 
differentiation with respect to the field V. 
We will calculate the expression 



d 



5T 



din A 8V y SV% 



(65) 



p=0 



We note that the regularization by higher covariant derivatives is essentially used here, 
because it allows differentiating the integrand and taking the limit of zero external mo- 
mentum. 

Using Eqs. lj5T|) . (JMJ), and |H2l) it is easy to see that 







d In A 5 V b y 



tr 



(T 



a , t mD 2 x f S 2 T 



IQd 2 \Scf> 



-i 



m „ mD 2 



8 2 T 



e 2 C(i?)^9 2 n 1/2 ^^ J 



IQd 2 \54>i84>i 
d 4 q 1 d f m 2 J \ 



-i 



+ h.c. 



p=0 



(66) 



(All similar expressions are written in the Euclidean space after the Weak rotation.) 

Let us proceed to the calculation of all other contributions in Eq. (|5fl|) . For this 
purpose we note that they contain the sum of subdiagrams, presented in Fig. [TJ The 
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y x z y x 

Figure 1: The sum of subdiagrams in the Schwinger-Dyson equations. 

external line in these subdiagrams corresponds to the background field V, and expressions 
for propagators were presented above. After simple transformations it is easy to see that 
expressions for the sum of these subdiagrams in the limit p — > are written as 



D 2 D 2 (v ° D2yD2y 5 s Uv GD2yL>2y 5 s 

y y \ y Q4d 2 (d 2 G 2 + m 2 J 2 ) yz J y A(d 2 G 2 + m 2 J 2 ) yz 

D 2 / \ CD 2 

d 8 x j^5 8 xy (2d 2 U 1/2 V x + iDj^DV^ 



-6 8 - 

d 2 G 2 + m 2 J 2 xz 



nn n 2 en 2 
-D a \ a z 5 8 - D 2 \ —5 5 8 (67) 

y 2(d 2 G 2 + m 2 J 2 ) yz y A{d 2 G 2 + m 2 J 2 ) yz 



or 



„ _ „ / mJD 2 \ mJD 2 _ 
= DDI V - S 8 I + V - S 8 = 

y A y lQd 2 (d 2 G 2 + m 2 J 2 ) yz J y d 2 G 2 + m 2 J 2 yz 

1 mJD av D 2 D 2 _ 1 , mJD 2 D 2 

8 ^ 2 (9 2 G 2 + m 2 J 2 ) V 16 y v & i (d 2 G 2 +m 2 J 2 ) yx 1 ; 

depending on the chirality of the ends. Terms in subdiagrams, containing the integral 
over d 8 x, have chiral projection operators, acting on the ends y and z. Due to identities 
(|60|) . this allows using simpler vertexes (JHTJ) and (|62|) for the calculations. After simple, 
but rather long calculations we obtain that the contribution of these terms to Eq. (|65|) is 



fdR) 6- #n 1/2 ^ / J^i A ln + m . j.) . (69) 

It is also easy to see that a contribution of the terms, which are proportional to D 2 V in 
Eqs. (|67|) and (|68|) . is 0. However, a contribution of the terms, proportional to D a Y is 
nontrivial and contains the unknown functions / and h. It can be written as 








«91nA 



tr / d 8 xd 8 yV y D a V C[ 



D n7 D 



az^ z 



5 3 T 



Aid 2 5]+5V 



z=x,p=0 



(70) 

D a V V 

where the symbol I a (we use the notation in [9] here) means that in the expression for 
the exact propagator it is necessary to make the substitution 
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D 2 D 2 5 8 XZ -> 2D 2 D a 5 8 xz . (71) 
After substituting propagators and vertex functions the considered contribution will be 



8 f d A q ( rn 2 j(2q 2 J'-(J-l)^ l 6G f + 16m 2 Jh" 

(72) 

Collecting all contributions to the effective action, we obtain the following result 
[Z = 1, a dependence on Z will be restored later): 



e C(R) 5 d U 1/2 5 xy — J — ^ - g 2 G2 + m2j2 J 



9 5 2 r 



01nA5V*5V£ 



+m 2 J 2 



^ = . . . + eV(H) ^H V2 ^ ^ / (in + 

m 2 J \ ™ 2 J (V J' ~ {J ~ 1)) 16G/ + 16m 2 J/i 



g 2 G 2 + m 2 J 2 J fq2Q2 + m 2 j2 j g 2 G 2 + m 2 J 2 

—similar terms for the Pauli-Villars fields >, (73) 



where dots denote contributions of the gauge fields and ghosts. Comparing the obtained 
result with the corresponding expression for the supersymmetric electrodynamics, we see 
that they differ in the multiplier C(R). Nevertheless, one more difference is possible: 
In the electrodynamics explicit calculations show that there is an identity, which can be 
written as equality to of contribution l|7T))l. The problem, if a similar identity is valid 
in the non-Abelian case, is not yet solved. In general, it is not necessary, because there 
are new diagrams, which contain a matter loop and vertexes with self-action of the gauge 
superfield, in the non-Abelian case. That is why in this paper we will not assume existence 
of such identity. 



4 Exact Gell-Mann-Low function. 



Investigating the supersymmetric electrodynamics |Hj reveals that there are diver- 
gences only in the one-loop approximation if the higher derivative regularization is used, 
while the Gell-Mann-Low function coincides with the exact NSVZ /3-function. Let us 
consider a non-Abelian theory. 

The Gell-Mann-Low function is defined by the dependence of finite part of the two- 
point Green function on the momentum in the limit m — > 0. And so, let us consider the 
massless case and write the renormalized two-point Green function as 

if = -y^tr J ^ V(-p)d 2 Tl 1/2 V(p)d-\a^/p), (74) 

where a is a renormalized coupling constant. It is known [213 EI] that the function d 
is invariant under renormgroup transformations, i.e. does not depend on a choice of 
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the renormalization prescription. The Gell-Mann-Low function, which we will denote by 
(3(a), is given by 



md(a,ii/p) 



d 



<9 hip 



d(a, \ijp). 



(75) 



Because the function d is renorminvariant, the function (3(a) will not also depend on 
a particular renormalization prescription. This is its difference from the function b(a), 
which is constructed from the renormalized coupling constant by the differentiation with 
respect to the logarithm of the normalization point. 

In order to obtain physical results it is necessary to calculate the function d in the 
massive case, and then to impose the boundary condition 



d(a, fi/p, mj \x) 



p=0 



(76) 



where a p h is an experimentally measured coupling constant, at a given fi. Moreover, it is 
necessary to impose one more boundary condition to the Green functions of the matter 
superfield: 



G(a, fi/p, m/ (i) 



m 



ph 



J(a, /i/p, raj /i) 



(77) 



"ph 



in which m p h denotes a physical mass. (This equality assures that the propagator of the 
matter superfield will have a pole at the physical mass.) The boundary conditions written 
above allow expressing the parameters a and m (at a given fi) through the physically 
observed values. Nevertheless, the results of this paper do not allow finding explicit form 
of the functions d, G, and J in the massive case. And so, we can not relate a and m with 
the experimentally observed values. However, it is possible to obtain an expression for 
the contribution of matter superfields to the Gell-Mann-Low function: 

If d denotes the function d calculated at Z = 1, the following equation takes place in 
accordance with Eq. ([731) 







■d 1 (a ,A/p) 



x 



d\np 

d 4 q 1 d 
(2tt) 4 V<¥ 
M?J f 







p=0 



<91nA 



d 1 (a ,^/p) 



-16tiC(R) 







p=0 



<91nA 



-E 



c;- 



ln(g 2 G 2 ) + X(g 2 /A 2 ) -$>m [q 2 G 2 PV + Mfj[ 

i 

-J2^X( q 2 /A 2 , ai )\, 



where the dimensionless functions X are obtained from the equations 



(78) 



<9X(g 2 /A 2 ^ 
dq 2 

dX(q 2 /A 2 ,a i 
dq 2 



16/ 



G 



M 2 j(2q 2 J'-(J-l)) 16Gf+16M 2 Jh 



q^q 2 G 2 + M 2 J 2 ) q 2 G 2 + M 2 J 2 



(79) 
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in which Mj = g^A. 

The obtained integral is reduced to the total derivative in the four-dimensional spher- 
ical coordinates, only the substitution at the low limit being different from [8j: 



«91nA 



Jpv(0, a. 



= ^ ^ { In G(0) 2 + X(0) - £ c, In (m* J Py (0, a, 
^CiX(0, a. 



1 / 01nG 1 
TT 1 H <91nA 2,9mA 



(80) 



g=0 



Here we took into account that the functions X(q 2 / A 2 , <2j) and Jpv{q 2 /A 2 , e^) in the limit 
g — > tended to some finite constants, because they were defined by convergent (even in 
the limit A — > oo), dimensionless integrals, which did not contain infrared divergences. 
Because both parts of Eq. (|80|) depend on a and A/p, it allows finding the expression 
for d,Q 1 up to an insignificant numerical constant 



do Vo, A/p) = --C(R) fin — - In G(a , A/p) - J-X(a , A/p)) + const. 

TT \ p 2 / 



(81) 



Making substitution (|3~K|) in Eq. ([KU| it is possible to restore the dependence of the 
effective action on the renormalization constant Z: 



d 1 (a,fi/p)=d 1 (a ,A/p) + 



C(R) 



TT 



In Z(a, A///). 



(82) 



Differentiating this expression with respect to the momentum p, we obtain the Gell-Mann- 
Low function: 



8 d - = -Wi + ^ + i«y ,83) 



d\np tt \ d\np 2d\np / 

Using Eq. (|79l) . we rewrite this contribution of the matter superfield as 

C{R) ( 16p 2 f 
— (1 - 7(a) - hm n —F7- 

TT \ P^O G 



(84) 



where 7(a) is the anomalous dimension of the matter superfield: 

7^d(a,/i/p)) = —^-lnZG(a,/i/p). (85) 

(Expression contributes to the Gell-Mann-Low function and should be compared 
with the numerator of Eq. (JIJ.) Let us again mention that the obtained result does not 
depend on a choice of renormalization scheme. 

We note that the contribution of the last term in Eq. (J81J) was always in all diagrams 
calculated up to now. Moreover, if the identity 



din A J - — " '* 4^2 gj+gy 



= 0, (86) 

z=x,p=0 
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is valid in the non-Abelian case, this will be always true. Nevertheless, in the non-Abelian 
case this statement requires an additional verification. That is why we for the generality 
write the correction, corresponding to a nontrivial contribution of the function /. 

We note that the function / can actually affect only diagrams in three and more loops. 
Therefore, if such contribution appeared calculating the scheme dependent function b(a), 
it could be always eliminated by a proper redefinition of the coupling constant. 

Finally we note that due to the similarity of the result with the electrodynamics case, 
if the contribution of the function / is absent, there are divergences only in the one-loop 
approximation. Details of the corresponding calculation are presented in Ref. 0. This 
completely agrees with the structure of the anomalies supermultiplet and is a particular 
feature of higher covariant derivative regularization using |22|. 

5 Conclusion 

In this paper a contribution of the matter superfields to the two-point Green function of 
the gauge field is calculated by the Schwinger-Dyson equations and Slavnov-Taylor identi- 
ties for the supersymmetric Yang-Mills theory regularized by higher covariant derivatives. 
This regularization is very essential, because the method of calculating, used here, is based 
on the differentiation with respect to the regularization parameter and taking the limit 
of zero momentum. 

In many respects the obtained results are similar to the corresponding results in the 
electrodynamics. However, it is necessary to mention some differences: 

1. The calculation was made using the background field method, because this method 
allows preserving the manifest gauge invariance of the effective action. Without its using it 
would be also necessary to calculate the renormalization constant for the gauge superfield 
that make calculations essentially more complicated. 

2. There is the multiplier C(R), defined by Eq. (|2j) in the non-Abelian case. The 
same multiplier is also present in the exact NSVZ /^-function. 

3. Deviations of the result from the exact /3-function, due to the contribution of the 
function /, which can not be determined from the Slavnov-Taylor identity, are possible 
hypothetically. Such contributions seem to disappear in the electrodynamics due to a 
new identity for the Green functions. However, there are diagrams, containing vertexes 
with self-action of the gauge field, in the non-Abelian case, which can break the addi- 
tional identity. Contributions of these diagrams were not yet calculated explicitly in the 
considered regularization. That is why we also present the result for the /3-function, if 
there is a nontrivial contribution of the function /. 

We note that according to the above discussion, it is possible to point two possible 
reasons why the three-loop calculations with the dimensional reduction differ from the 
NSVZ /3-function: 

1. There is a nontrivial contribution of the function /. 

2. In preceding papers instead of the Gell-Mann-Low function, the /3-function, defined 
as a derivative of the renormalized coupling constant in the MS-scheme, is calculated. 

At present we can not affirm which of this variants (or both ones simultaneously) is 
realized. This requires either explicit calculating or finding an origin of identity (|86[). 
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Moreover, it is necessary to calculate a contribution of the gauge fields and ghosts 
using Schwinger-Dyson equations and Slavnov-Taylor identities. Now this work is in 
progress. 
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